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ABSTRACT 

We describe a new software package capable of predicting the spectra of solar-system 
planets, exoplanets, brown dwarfs and cool stars. The Versatile Software for Transfer of 
Atmospheric Radiation (VSTAR) code combines a line-by-line approach to molecular 
and atomic absorption with a full multiple scattering treatment of radiative transfer. 
VSTAR is a modular system incorporating an ionization and chemical equilibrium 
model, a comprehensive treatment of spectral line absorption using a database of more 
than 2.9 billion spectral lines, a scattering package and a radiative transfer module. 
We test the methods by comparison with other models and benchmark calculations. 
We present examples of the use of VSTAR to model the spectra of terrestrial and 
giant planet in our own solar system, brown dwarfs and cool stars. 

Key words: radiative transfer - techniques: spectroscopic - planets and satellites: 
atmospheres - stars: atmospheres - brown dwarfs. 



1 INTRODUCTION 

' Until recen tly the modelling of the atmospheres of stars (e.g. 
; iGravl bOOSl ) and the modelling of the at mospheres of the 
■ Earth and other solar-system planets (e.g. lLiou|[20o3 ) have 
' developed largely independently. Models of stars applied to 
, high temperature objects with effective temperatures Tcff > 
' 3000K, with opacity dominated by the line and continuum 
absorption of atoms and atomic ions, whereas planetary at- 
mosphere models applied to cool objects Tes ~ 100-300K 
where the important processes were molecular absoprtion 
and scattering from molecules and cloud particles. 

This situation changed with the discovery in the mid- 
1990s of the first un ambiguo us brown dwarf, G1229B. 
iNakaiima et al.l 1 19951 : Qppcnh eimer et al.l Il995l ) and the 
first hot Jupiter planets beginning with 51 Peg b 
l|Mavor fc Quelodl 19951 : iMarcv et al ]| 19971 ). Many more such 
objects have now been discovered and reveal that planets 
and brown dwarfs populate an intermediate range of temper- 
atures not explored previously. This has led to the require- 
ment to develop new methods to model these atmospheres 
that cover the effective temperature range from below lOOOK 
to more than 2000K. 

One widely used approach has been to adapt stel- 
lar atmosphere codes to handle the lower temperatures 
encountered in exoplanets and brown dwarfs. Models of 
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this type are descr ib ed for example bv iTsuii et al.l 
lAllard et all (|200ll). [Barm an, Hauschildt fc AUardl 
and iBurrows. Sudarskv fc Hubenv (2003, 2006). An alter- 
native approach, and the one we follow in this paper, 
is to take models originally used for the atmospheres of 
the Earth or other solar system planets and adapt them 
to handle the higher temperatures needed for exoplan- 
ets and br own dwarfs. S uch an ap proach has been de - 
scribed by iMarlev et all l|2002l ) and jPortnev et al] (|2005l ) 
who use a model based on one origin ally used to mode l 
the atmospheres of Tit an and Uranus (jMcKav et al ] |l989l : 
iMarlev fc McKavll 19991 ). 

One important difference between the techniques used 
in stellar atmosphere modelling and in Earth atmosphere 
modelling is the approach to radiative transfer. In stellar 
atmospheres simplified treatments of scattering, such as the 
assumption of isotropic scattering, are commonly adopted. 
This is justified by the fact that scattering is in most cases 
not an important source of opacity in stellar atmospheres, 
and where it does become significant, in the form of Rayleigh 
scattering from molecules in cool stars, and scattering from 
electrons in hot stars, the phase functions are forward- 
backward symmetric. Scattering from clouds and aerosols in 
the Earth atmosphere can, however, result in strongly for- 
ward peaked phase functions, and to properly model such 
cases radiative transfer methods that more rigorously han- 
dle multiple scattering with anisotropic phase functions are 
needed. While these techniques build on the classic work of 
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astronomers such as Chandrasekhar and van de Hulst, much 
recent development of such method s has been in the context 
of Earth atmosphere research (e.g. lLioul l2002''l. 

Clouds are now known to be important not just in 
the atmospheres of all the solar system planets, but in 
many brown dwarfs and even in late M dwarfs. Thus 
for all these objects a more rigorous treatment of radia- 
tive transfer is desirable. Such an approach is particu- 
larly important for modelling the reflected light from ex- 
oplanets, because the angular dependence of scattering is 
an important factor in deter mining the phase variation 
around a planet's orbit (e.g. Seager. Whitnev fc Sasselovl 
l2000l : ICahov. Marlev fc Fortnevll20ld~ 

In this paper we describe the methods used in the VS- 
TAR (Versatile Sofwtware for Transfer of Atmospheric Ra- 
diaition) atmospheric modelling software. VSTAR was orig- 
inally developed as a way of modelling the spectra of the 
sol ar system pla nets, and an early version of it is described 
bv lBailevI l|2006l ). In this paper we describe the current ver- 
sion of VSTAR which can now handle a wide range of at- 
mospheres ranging from those of the coolest solar system 
planets up to stars with temperatures of '^SOOOK, and thus 
including the brown dwarfs and hot Jupiter type exoplanets. 



2 THE VSTAR MODEL 

2.1 Ionization and Chemical Equilibrium 

The Ionization and Chemical Equilibrium (ICE) package of 
VSTAR is used to determine the equilibrium chemical com- 
position of an atmospheric layer given its elemental abun- 
dances, pressure and temperature. ICE handles gas phase 
chemistry, ionized species and the formation of solid and 
liquid condensates. Full details of the methods employed 
in this package will be described elsewhere so only a brief 
description is given here. The te c hniques used are si milar 
to those described bvlTsuiil (|l973h . lAllard et al.l (|200ll ) and 
iLodders fc FeglevI l|2002l ). For each compound considered in 
the model the equilibrium constant of formation Kf from the 
elements is required. These are used in a set of equations for 
the mass balance of each element, and for the charge bal- 
ance to solve for the abundances of each molecular species. 
Kf is related to the Gibbs free energy of formation A/G° 
through: 



AfG° = -RTlnK 



f 



(1) 



Where R is the gas constant. This relationship is in- 
dicative of the link between the technique we use based 
on equilibrium coinstants, and the alternative technique 
for chemical models based on m inimization of the total 
Gibbs free energy of t he system (IS harp fc Huebneil Il990l : 
ISharp fc Burrows|[2007l ') Both AfG° and Kf are functions of 
temperature (T) and are available in standdard compilations 
of thermochemical data. Our main source of thermochemi- 
cal data was the fourth edition of the NIST-JANAF therm o- 
chemi cal t ables (Ch ase 1998) . We used data from iLodders 
l|l999h and lLodderl (|2004l ') for several compounds for which 
the JANAF data have been shown to contai n errors. Ad - 
ditional sources of thermochemical data were iBarinI jWd^ ) 
and lRobie fc Hemingwavl (119951 ). 



Some important species, however, do not have data in 
any of these tabulations, and for some others the available 
data do not extend to sufBciently high temperatures for stel- 
lar atmosphere models. For a number of gas-phase species 
we have therefore calculated our own thermochemical data. 
To do this we first calculate the partition function using the 
spectroscopic constants of the molecule. We use a rigid rota- 
tor harmonic oscillator model (with additonal corrections for 
anharmonicity and centrifugal distortion) for the rotational 
and vibration al levels, and direct su mmation over the elec- 
tronic levels (|Maver fc Mavej Il940l ). The thermodynamic 
functions can then be derived from the partition function 
and the dissociation energy of the molecule. The techniques 
are essentially the same as those use d in construc ting the 
NIST-JANAF thermochemical tables l|Chaselll998l ') for gas 
phase species, and we have tested our methods by repro- 
ducing results published in these tables. Species we have 
calculated new thermochemical data for include FeH, CrH, 
CaH, TiH, RbCl and RbF. 

Our thermochemical model currently includes 143 gas- 
phase and condensed-phase compounds of 27 elements. 
While this is not as extensive as some other models, we 
have been careful to include all species listed a s imp ortant 
in previous studies such as iBurrows fc Shard l|l999t ). Our 
chemical model predictions s how good agreement with pre- 
vious results such as those of ILodders fc FeglevI (12002') and 
ISharp fc Burrow s (2007). Figure[l]shows the results of calcu- 
lations of the mixing ratios, or mole fractions, for a number 
of important species at 1 atmosphere pressure over a range of 
temperatures. It can be directly compared with figure 17 of 
ISharp fc Burrowd \200f). The model can easily be expanded 
by adding new species to our thermochemical database. 

Our current model assumes true equilibrium chemistry 
in that when condensates form, the condensed phase is as- 
sumed to remain in equilibrium with the gas phase. This 
assumption is commonly made in such models. However, 
it may not be true in real cases since condensed material 
can fall under gravity and rain out of the system. Models 
that take a ccount of this "rainout" proces s are described by 
iMarlev et al.. (,200^) and iFreedman et al.l ([2008.1 



2.2 Molecular Absorption Lines 

Absorption lines due to rovibrational and electronic transi- 
tions of molecules are the most important features of the 
spectra of planets, brown dwarfs and the coolest stars. For 
use with VSTAR we have collected a line database that cur- 
rently contains more than 2.9 billion lines. The line lists used 
with VSTAR are listed in table[T] The line parameters are in 
a variety of different formats. The essential data needed for 
all lines are the line position (as wavelength or wavenum- 
ber), the line intensity, and the lower state energy of the 
transition (ususally given in cm~^ or eV). 

Line intensity is usually given in one of three forms. 
HITRAN, GEISA (see section 12.2. ip and similar lists give 
a line intensity So in units of cm molecule"^ (sometimes 
given as cm'^ mol~^ cm~^) at a reference temperature To- 
For HITRAN To is 296 K, but some other lists use a different 
reference temperature. 

T he line intensity at te mperature T is then calculated 
using (|Rothman et al lliggsl ): 
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2.2.1 HITRAN and GEISA 
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Figure 1. Mixing ratios on several important species as a func- 
tion of temperature for a solar composition gas at 1 atm pressure 
calculated with our chemical equilibrium model. The r esults can 
be compared with figure 17 of lSharp fc Burrows! |20o3) 



_ SoQ{To) exp{~C2Ei/T) [1 - exp(-e2i/o/r)] 

Q{T) exp{-C2Ei/To) [1 - exp(-C2//o/ro)] ^ ' 

where vo is the line frequency in cm~^, Ei is the lower 
state energy in cm~^, C2 is the second radiation constant (= 
hc/k) and Q{T) is the partition function (or total internal 
partition sum). 

Line intensities may also be quoted in the form of Ein- 
stein A coefficients (^21 )• The line intensity S (in cm/mol) 
at a temperature T can be calculated from A21 using 
ISimeckova et al.ii2006iy 



S = 



5^2 



exp(-C2£i/T)(l-exp(-C2i'o/r)) (3) 



Here g is the statistical weight of the upper level of the 
transition. In most cases g can be written as (2J + l)gs 
where J is the rotational quantum number of the upper 
state, and gs is the nuclear spin degeneracy. Different for- 
mulae are n eeded if hyperfine stru cture is included in the 
line list (see ISimeckova et aLlfiooil ). The gs used in equa- 
tion [3] needs to be consistent with that used in calculating 
the partition function (see section [!^. 5. 2 1) . 

In astronomical line lists, such as those of iKurucj 
(|2005l ). the oscillator strength /, usually tabulated as gf or 
log{gf) (where g is the statistical weight) is the usual form 
for listing line strengths. The line intensity at temperature 
T can be calculated from g f using; 



HITRAN (High resolution TRANsmission) and GEISA 
(Gestion et Etude des Informations Spectroscopiques Atmo- 
spheriques) are compilations of molecular spectroscopic pa- 
rameters designed primarily for the Earth atmosphere, but 
often useful for the atmospheres of other planets. Both are 
updated every few years with the most recent releases be- 
i ng HITRAN 2008 (iRothman et al.ll2009l ) and GEISA 2009 
( Jac guinet-Husson et al.ll2008h . The line parameters are gen- 
erally best suited to low temperature models and more com- 
plete line lists are usually needed for higher temperature at- 
mospheres such as those encountered in bro wn dwarfs, hot 
Jupiter type planets and stars. For example iBailevI (|2009| ) 
showed that HITRAN and GEISA were significantly incom- 
plete for water vapour line parameters at temperatures of 
500- TOOK encountered in the Venus lower atmosphere. 

2.2.2 HITEMP 

The HITEMP database is a companion to HITRAN con- 
taining line data suitable for use at higher tempera- 
tures. Th e latest editi on of HITEMP was released in 2010 
( Roth man et al.l l2010ll and repla ces an earlier edition de- 
scribed bylRothmaneFal] l|l995h . HITEMP contains data 
for five species, H2O, CO2, CO, NO and OH. The data for- 
mat is the same as that used in HITRAN, and HITEMP is 
consistent with HITRAN, in the sense that lines common 
to both databases have the same line parameter values, al- 
though HITEMP includes many lines not in HITRAN. 



2.2.3 Water Vapour — H2O 

Water vapour is an important absorber in atmospheres 
ranging from stars to the terrestrial planets. A num- 
ber of line lists are available for water vapour at 
high temperatures. These include the HITEMP list 
(Rot hman et all I2qi0h. the SCAN list djorgensen et all 
I2OO1I) . th e list of [Partridge fc Schwenkd (|l997h and the 
BT2 list ijBarber et all 120061). Compari s ons o f the various 
lists h ave b een reported b y Allard et al.l (j2000l ) , Ijones et al.l 
(!2003f) and iBailevI (|2009l '). We normally use the BT2 list 
which is the most exte nsive and is based on the ac urate 
dipole moment surface o f ISchwenke fc Partridge! (j2000l ). The 
HITEMP list for H2O is based on BT2. 

BT2 includes only th e main isotop o logue of water 
(H2^'*0). The VTT list of IVoronin eTaP (|2010l l provides 
a similar list for HDO. In the atmosphere of Venus deu- 
terium is enhanced over terrestrial abundances by a fac- 
tor of 100-150 so HDO absorption is significant. Line data 
for a range of other isotopologues of water at 296K and 
lOOOK are available from the Spectra Information Sys- 
tem jhttp:/ /spectra iao.ru[ ) at the Institute of Atmospheric 
Optics, Toms k. These data are based on the analysis of 
iSchwenke fc Partridge! (|2000l l. 



where e is the electron charge, is the electron mass 
and other symbols are as above. 



2.2.4 Carbon Dioxide — CO2 

Carbon dioxide is an important constituent of terrestrial 
planet atmospheres. Until recently it has not been thought 
to be important in giant planet atmospheres. However, re- 
cent results from transitting extrasolar planets (|Swain et al.l 
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Table 1. Molecular Line Lists used with VSTAR 



Molecule 


List 


Number of Lines 


Reference 




39 molecules 


HITRAN 2008 


2,713,968 


Rothman et al. (2009) 




50 molecules 


GEISA 2009 


3,807,997 


Jacquinct-Husson et al. ^008) 




H20 


BT2 


505,806,202 


Barber et al. (2006) 




H20 


SCAN 


101,455,143 


j0rB:ensen et al. (2001) 




H20 




65,912,356 


Partridge & Schwcnke (1997) 




H20 


HITEMP 


111,377,777 


Rothman et al. ("2010") 




HDO 


VTT 


697,454,528 


Voronin et al. (2010) 




CO2 




7,088,178 


Pollack et al. ( 1993) 




CO2 


HITEMP 


11,377,777 


Rothman et al. (2010') 




CO2 


CDSD-296 


419,610 


Tashkun et al. (2003) 




CO2 


CDSD-1000 


3,950,553 


Tashkun et al. (2003) 




CO2 


CDSD-Venus 


11,730,277 


Tashkun et al. (2003) 




CO 


HITEMP 


115,218 


Rothman et al. (2010) 




CO 




134,421 


Goorvitch (1994) 




CH4 (cool) 




339,690 


see text (section |2. 2. 6|| 




CH4 (hot) 




134,862,336 


see text (sectionl2. 2.611 




NH3 




3,249,988 


Yurchenko et al. (2009) 




NH3 


BYTe 


1,138,323,351 


Yurchenko, Barber & Tennyson 


(2011) 


TiO 


SCAN 


12,837,150 


Jargensen (1994) 




TiO 




37,744,499 


Schwenke (1998) 




TiO 




11,369,552 


Plez (1998) 




VO 




3,171,552 


Plez, B., private communication 




CaH 




124,615 


see text (section 12. 2. 9[ 




MgH 




23,315 


Week et al. (2003): Skorv et al. f 


2003) 


MgH 




119,167 


Kurucz (2005) 




FeH (F-X) 




116,300 


Dulick et al. (2003) 




FeH (E-A) 




6,357 


Harereaves et al. (2010^ 




CrH 




14,255 


Burrows et al. (2002) 




TiH 




199,073 


Burrows et al. (2005) 




CH 


SCAN 


114,567 


^grgsnsen et al. (1996) 




CN 


SCAN 


2,245,378 


j0rBensen & Larsson Q990) 




C2 




360,887 


Ouerci, Ouerci & Kunde (1971): 


Ouerci. Ouerci & Tsuii (1974) 


HCN/HNC 




34,433,190 


Harris et al. (2006) 



l2009al lb[) suggest CO2 may be important in these atmo- 
spheres as well. Recently CO2 has also been det ected 
in brown dwarfs (|Yamamura. Tsuii fc Tanabd I2OI0I ). At 
low temperatures HITRAN, GEIS A or CDSD-296 (Ca rbon 
Dioxide Spectroscopic Databank, iTashkun et al.|[2003l ) can 
be used. H igh temperature lin e lists include the list de- 
scribed by iPollack et all (| 19931 ) and the CDSD-1000 list 
l|Tashkun et all |2003| )." The latter list contains more than 
3 million lines and is complete to an intensity of 10"'^^ cm 
mol"^ at 1000 K. 

The dense CO2 atmosphere of Venus provides a par- 
ticular challenge to the modelling of CO2 absorption with 
pressures up to 90 bars of almost pure CO2 and requires a 
line list with a much d eeper intensity cutoff than is normally 
necessary. The list of iPoUack et al.l (|l993t ) is the one that 
has normally been used for modelling the Venus deep atmo- 
sphere. However, an alternative is the CDSD-Venus list, a 
version of the CDSD list with a deeper cutoff of 10~^° cm 
mor^ at 750 K. 

The HITEMP (jPothman et al.l l2010l ) list for CO2 is 
based on the CDSD line lists. 



2.2.5 Carbon Monoxide — CO 

Carbon monoxide is significant in terrestrial planet atmo- 
spheres, and is also observed in Titan. It also becomes 



important in high temperature atmospheres where equilib- 
rium chemistry favours CO over CH4. For high temperatures 
the HITEMP lin e list (|Rothman et al.ll201(i ) or the list of 
iGoorvitchI (|l994l ) can be used. 



2.2.6 Methane — CH4, 

Methane is an important absorber in the atmospheres of the 
solar system giant planets and in Titan. It has b een detected 
in th e atmospheres of extrasolar giant panels (|Swain et al.l 
l2008h . It is also important in brown dwarfs, the presence of 
methane absorption features in the near-IR being the defin- 
ing characteristic of the T-dwarf class. 

However, methane has a complex spectrum due to the 
presence of coincidences between its four vibrational modes 
that result in a series of interacting states known as polyads, 
spaced at intervals of about 1500 cm"^. Most line data 
is either based on experimental measurements or on effec- 
tive Hamiltonian models, which exist only for the lowest 
few pol yads. The Spherical Top Da ta System (STDS) soft- 
ware o f lWenger fc ChampionI (| 19981 ) can be used to generate 
line lists from these effective Hamil tonian models. A g ood 
model has recently been obtained (|Albert et al.l |2009| ) for 
the ground-state and the lowest three polyads (the Dyad, 
Pentad and Octad) and this allows reliable prediction of the 
low temperature methane spectrum from - 4800 cm~^. 
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This model is used in th e 2008 update of the m ethane line 
parameters in HITRAN l|Rothman et al.llioogl 'l. A prelimi- 
nary model is avai lable for t he next polyad, the Tetradecad 
l|Boudon Rev fc Lo otc 200^; [Robert et ai]|200ll ). 

Above 4800 cm~^ many of the line parameters included 
in HITRAN are emp irical measure ments at room tempera- 
ture originating from lBrownl ()2005l ). These mostly lack lower 
state energies (HITRAN lists fictitious lower state energies 
of 555.5555 or 333.3333) and so cannot be used to derive re- 
liable line intensities at other temperatures. However, some 
lines in the 5500 - 6150 cm~^ region have empirical lower 
state energies derived from me asure ments at multipl e tem- 
peratures from iMareohd (|l990l ) and lGao et all (|2009[ ). 

Recently much improved data for the low tempera- 
ture methane spectrum over the range 1.26 - 1.71 
(5852 — 7919 cm~^ has become available from labora- 
tory m e asurements d e scribe d bylWang, Kassi fc Campargue 
tOld): IWang et all (l201lh . ICampargue et al.i |2O10l l and 
iMondelain et al.l (|201ll ')~These data are based on a com- 
bination of room temperature and cryogenic (^80K) mea- 
surements that allow empirical determinations of lower state 
energies, and includes deep cavity ring down spectroscopy 
of the weak lines in the window regions between the strong 
absorption band systems. An ap plication of this l ist to the 
spectrum of Titan is reported bv lde Bergh et alT] (|201ll ). 

For low temperature atmospheres, such as the solar sys- 
tem giant planets or Titan, we use a composite list that com- 
bines the la boratory data described a bove, as well as data 
described bv lNikitin et al.1 (|201G| . |2011| 'I. with HITRAN data 
from - 4800 cm"\ We use STDS to fiU in the remaining 
gap in coverage between 4800 and 5500 cm"'^. An early ver- 
sion of this l i st and its application to Titan is described by 
iBailev et al.l (1201 ll ) . The current version of the l ist is further 
impro ved by including the most recent data of I Wang et al.l 
l|201ll) andlMondclain et al. (2011). 

For high temperature atmospheres, such as those of 
brown dwarfs and hot Ju piters, we use a line list com puted 
with the STDS software (jWenger fc Championiri998l ) up to 
J=60 for most band systems, and up to J=50 for the most 
complex system modelled (the Tetradecad-Pentad). This 
gives a list of nearly 135 million lines . Beca use the effec- 
tive Hamiltonian model of I Albert etall l|2009l ) involves high 
order polynomial fits to empirical data usually limited to 
the range up to J=20-30, it is not suitable for extrapola- 
tion to these high J values. We inst ead use the low orde r 
effective Hamiltonian parameters of iBorvsov et al.l l|2002l ) . 
These provide a poorer fit to observed line positions ('^ 0.5 
cm~^) but allow more reliable extrapolation to high J. For 
most purposes we use a smaller list extracted from the full 
list with a line intensity cutoff of 10~^^ cm moP^ at 1500 
K giving a list of 15.8 million lines. 

While this approach to the hot methane spectrum 
is similar to that used by other g roups modelhng brown 
dwarf and exoplanet spectra (e.g. [ Freedman et al.l I 20081 : 
iHomeier et al. I I2OO3I : iBorvsov et al.ll2002l ) it is important to 
understand that while these models include high rotational 
levels, they do not include higher vibrational levels that are 
needed at these temperatures and many hot band systems 
are therefore missing. Methane band strengths will therefore 
be underestimated by amounts that will increase for shorter 
wavelengths and higher temperatures. 

Above 6500 cm~^ there is little useful methane line data 



for high temperatures. The low temperature line list can be 
used but will omit both hot bands and high J transitions. 
An alternative is available in the form of low resolution 
absorption coefficients , band models or k-distribution pa- 
rameters for methane ((strong et al.lll99j : llrwin et al. I I2OO6I : 
iKarkoschka fc Tomaskoll2010l ). The first two sets cover the 
near-IR, and the latter extends through the visble as well. 
However, these data are not in an ideal form for use with 
line-by-line programs such as VSTAR and are not designed 
for high temperature use. The optimal way to use such 
datas ets is in conjunction with the correlated-k method 
(Goodv et al.|[r98ai 



2.2.7 Ammonia 



NH3 



Ammonia is present in the atmospheres of Jupiter and Sat- 
urn, and is also found in cool brown dwarfs where it is seen 
most easily through a feature at 10.5 fim (jCushing et al.l 
I2OO6I). The spectrum of ammonia, has until recently, pre- 
sented similar problems to methane. The line parameters in 
HITRAN only include lines up to 5295 cm"'^ and have not 
been updated since the 2000 edition. However, recently a 
new line list (the BYT e list) for hot NH3 (up to ~1500 K) 
has been published bv lYurchenko et al.l (|201ll ). It contains 
more than 1.1 billion lines covering the range - 12,000 
cm~^. An NH3 list for use up to 300K and cove ring - 8000 
cm~^ was described bv lYurchenko et al.1 l|2009D . 



2.2.8 Metal Oxides — TiO and VO 

Bands of Titanium Oxide (TiO) and Vanadium Oxide (VO) 
are distinctive features of the spectrum of M- type stars. 
There are several l ine lis t available f or TiO from |j0rgensenl 
1 1994 ). ISchwenkd l| 19981 ) and iPlej (|l998l ). The lists cover 
five isotopologues o f the TiO mo l ecule. Comparison s of the 
lists can be found in lAllard et al.1 (|2000l ) and Pavle nko et al.l 
(|2006l ). The VSTAR database also includes a VO hue hst 
provided by Plez, B. (private communication), which is part 
of the line lists assembled by iGustafsson et al. I iliooai . 



2.2.9 Metal Hydrides — CaH, MgH, FeH, CrH and TiH 

Bands due to electronic transitions of metal hydrides are 
seen in M-type stars a nd some cooler ob jects . A line list 
for M gH is described bv lWeck etall (|2003h and lSkorv et al.1 
(l2003h . This list, and a similarly formatted list for CaH 
(which does not appear to have a published description) 
can be found on the web site of the University of Georgia 
Molecular Opacity ProjeclQ. An alternative MgH list which 
includes two a dditional MgH isotopologues is available from 
iKurucd (|2005l ). 

L ines lists are also ava ilable for Fe H llDulick et al] 



CrH (|Burrows et al.ll2002l ) and TiH l|Burrows et al.l 



2003 ). 



20051 ) 



These line lists are all similarly formatted. They cover only 
the most abundant isotopologue of each molcule, though 
methods for determining line position s for other isotop- 
logues are described. The iDulick et all (|2003l ) line list for 
FeH covers the F-X band, but absorptions in the E-A band 



http: / / www. physast . uga. edu / ugamop / index, html 
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Figure 2. Monochromatic absorption (in cm^ molecule ^) of 
CaH at a temperature of 1800 K and a pressure of 10 atm 



are also important in the 1.6 fim region. These are avail- 
able from a recent emp irically based line list described by 
iHar greaves et all (|201oh . 



2.2.10 CH, CN and Ca, HCN/HNC 

These bands are present in cool stars and become par- 
ticularly strong in car bon st ars. Line lists availabl e from 
|j0rgensen et al.l (Il996l) ( CH) [j^ rgensen fc LarssonI (|l990l ') 
(CNl.lQuerci et al l l| 1971] ) and lQuerci et all l| 19741 ) (C^ ) and 
iHarris et all (|2006l ) (HCN/HNC) are included in the VS- 
TAR database. 



2.2.11 Absorption Plots 

ISharp fc Burrows! |20o3) have presented plots of the 
monochromatic absorption (in cm^ molecule"^) for many 
molecular species present in brown dwarf or giant exoplanet 
atmospheres. We have made similar plots for the molecu- 
lar species described above. In all but one case our plots 
look e ssentially the same as those given bv lSharp fc Burrows! 

The exception is CaH where our absorption plot, 
given in figure [21 shows the well known absorption bands at 
0.63 8 and 0.683 ^m in M dw arfs, which are absent in figure 
9 of !Sharp fc Burrows! l|2007l ) for the same temperature and 
pressure. 



2.3 Atomic Lines 

While molecular lines are most important in the coolest ob- 
jects, lines of the alkali metals are significant in brown dwarfs 
and exoplanets, and many other atomic lines become impor- 
tant in M stars. VSTAR can make use of atomic line data 



in the format of the Kurucz line list s, or that of the Vi- 
enna Atomic Line Database (VALD, !Piskunov et al.![l995! : 
iKupkaet aDll999! ). 



2.4 Line Profiles 

The default treatment of spectral line profiles in VSTAR is 
to use a Voigt line shape in t he line cor e , app roxi mated us- 
ing the methods described bv lHumlicekl l| 19821 ) and lSchreieJ 
( 1992), out to 20 Doppler widths from the line centre. In the 
line wings a van-Vleck Weisskopf pr ofile is used which has a 
line s hape function 4'{v) of the form l|van Vleck fc Weisskopll 
1 19451 ): 



+ 



7 



(5) 



Where v is frequency, vo is the line centre frequency, 
and 7 is the line half width. In most cases this is essentially 
identical to a Lorentzian, but this profile provides the asym- 
metry that is needed at microwave wavelengths or with very 
broad lines, and goes to zero at zero frequency. 

In general it is found that continuing this profile out 
to large distances from the line centre overestimates the ab- 
sorption in the far wings. In cases where the precise far wing 
behaviour is important, VSTAR allows the line wings to 
be modified by multiplying by a correction factor, usually 
called a x factor. For example the far line wings of CO2 
are known to be sub-Lorentzian JPerrin fc Hartmannl!l98^ : 
iTonkov et al. |[l996,), and the correct profile is important in 
modelling the spec tra of the "windows" in t he Venus night- 
side spectrum (e.g. !Meadows fc Crisd[l99^ ). 

Where the precise line wing shape is not critical, it is 
usually sufficient to cut off the line profiles at a distance from 
the line centre, which can be chosen according to the typical 
line strengths and widths but usually ranges from 20 to 100 
cm~^. This approach improves the speed of computation. In 
either this case, or the x factor case, the line intensities are 
adjusted to ensure that the integrated line strength is not 
changed from its original value. 

The line profile requires a value for the pressure broad- 
ened hne half width (7 in equation [S]). HITRAN, GEISA 
and similarly formatted files such as HITEMP and CDSD 
include line width data in the form of an air broadened line 
half width coefficient 70 (in cm~^ atm~^, quoted for the 
reference temperature To), a self broadened line half width, 
and a temperature exponent n for the line width such that: 



l{T) = 70 



(6) 



Line widths for different transitions of a molecule vary, 
and are usually found to depend primarily on the rotational 
quantum numbers of the levels involved in the transition. 
Variation with the vibrational quantum numbers is much 
smaller, and is often ignored in empirical models for the line 
shape. 

Line widths depend on the type of gas responsible for 
the broadening. Hence the air broadened widths given in 
HITRAN type files are not usually what is needed, as we 
generally want the line widths for broadening in H2 and 
He (for giant planet and brown dwarf atmospheres) or CO2 
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Table 2. Relative broadening effects of different broadening gases 



Lines 
of 


H2 


He 


N2 


Broadened by 
O2 CO2 


Air 


H2/He 


References 


NH3 


0.832 


0.370 


1.0 


0.693 




0.939 


0.818 


Brown & Peterson (1994) 


CH4 


1.017 


0.643 


1.0 


0.943 




0.989 


0.961 


Pine & Gabard f2003) 


CO 


0.85 


0.64 


1.0 


0.90 


1.30 


0.98 


0.82 


Burch et al. (1962), Hartmann et al. (19881 


H2O 


0.777 


0.221 


1.0 


0.661 


1.50 


0.932 


0.694 


Gamacfie et al. fl996) 


CO2 


1.41 


0.59 


1.0 


0.81 


1.20 


0.962 


1.287 


Burch et al. ( 1969) 



(for terrestrial planets) . Unfortunately the data available on 
these alternate broadening gases is rather incomplete. 

Table [2] shows the relative line widths due to broaden- 
ing in different broadening gases, averaged over many transi- 
tions. The line widths are relative to that due to broadening 
in N2. Some general trends can be identified in the table. 
For example, helium always produces the smallest broaden- 
ing. O2 broadening is always less than that due to N2. CO2 
produces the largest broadening, However, H2 is sometimes 
less than, and sometimes more than N2. These trends can 
be used as a rough guide, when no broadening information 
is available for a specific broadening gas. 

Broadening for a mixture of gases can be obtained by 
averaging the broadening coefficients weighted by their par- 
tial pressures (or volume mixing ratios). Table[2]gives values 
for air (taken as 80% N2, 20% O2) and an H[2/He mixture 
typical of a giant planet or brown dwarf atmosphere (85% 
H2, 15% He) obtained from the other values in the table. 

VSTAR includes specific models for the line widths of 
several of the more important gases. For H2O broadened in 
CO2 we use the model described by Bailey ( 2009 ) based on 
the data of [Pelave et al.l (|198^ ) . For H2Q broa d ened in H2 
and He we use the data from iGamaclie et al.l l| 19961 ). For 
N H3 broadened by H2, He , N2 and O2 we use the results 
of lBrown fc Petersonl (|l994 ) for the dependence of width on 
the J and K q uantum numbers, and the temperature expo- 
nents are from iNouri et al.l l|2004 ) . For CO we use the model 
for the J depend ence from HITRAN (iRothinan et al.ll2005l ) 
for air, and from iRegaha-Jarlot et all (|2005l ) for H^;.. When 
CH4 line data are taken from HITRAN, the air broadened 
width from HITRAN is used for H2/He broadening also, 
since table [5] shows little difference in the two cases. For 
other sour ces o f CH4 line data the broadening is based on 
data from lPind j 19921 ). 

Atomic line lists usually give width data in the form of 
three coefficients for natural (7n), van der Waals (7™) and 
Stark (7a) broadening. Stark broadening is usually negligi- 
ble at the temperatures considered here. The line width 7 
needed for equation [5] can be calculated from; 



1 

47rc 



■yw{NH + C1NH2 + C2Nno) 



T 



10, 000 



+ 7n 



(7) 

where Nh, and TVhc are the number densities of 

atomic hydrogen, molecular hydrogen and helium. Ci and 
C2 are coefficients for the relative broadening effects of H2 
and He as compared with H, which can be derived from 
the polarizabilities of the different species. If 7™ and 7„ are 
no t available they are calc ulated using methods described 
bv lKurucz fc Avretj (|l98ll ). 



2.5 Partition Functions 

The calculation of line intensities from data tabulated in 
the line lists using equations [2] [3] or [4] requires the partition 
function Q{T) of the molecule or atom defined by. 



Q(r) = ^<;.exp(-|i 



kT 



(8) 



where Ei is the energy of level i relative to the ground 
state, gi is its statistical weight, and the sum is to be taken 
over all rotational, vibrational and electronic levels. 



2.5.1 Molecular Partition Functions 

VSTAR uses molecular partition functions from a range of 
sources. Polynomial approximations for the partition func- 
tions of many molecules for temp eratures from 1 000 K up to 
9000 K or h igher are given bv iSauval fc TatumI (1984) and 
llrwini (|l98ll ). A subroutine for calculating the partition func- 
tions of all molecules and isotopologues in HI TRAN is pro- 
vided with that database and is described by iFischer et al.l 
(2003^. The HITRAN partition functions are valid over the 
temperature range from 70 to 3000 K. An updated version 
of the H ITRAN partition fu nctions has recently been pub- 
lished bv lLaraia et al.l(|201ll ). but the work described in this 
paper used the earlier versions from iFischer et all l|2003h . 
Partition functions can be reconstruct ed from dat a in the 
NIST-JANAF thermochemical tables (|Chaselll99^ ). These 
tables do not list the partition functions explicitly, but the 
thermochemical properties listed for gas phase substances 
are derived from partition functions. The partition function 
can be calcu lated f rom t he tabulated data using methods 
described bv llrwini (|l988l ). 

Partition functions can also be calculated directly from 
spectroscopic constants using methods already described in 
section [27l] In some cases it is possible to calculate molecular 
partition functions by direct summation over energy levels as 
in the partition funct ions for H2O f r om[Vidler fc TennvsorJ 
(|2000h and CH4 from lWenger et al.l i20m . 



2.5.2 Partition Functions and Nuclear Spin 

In comparing the molecular partition functions from differ- 
ent source s it was noticed tha t the partition functions from 
HITRAN (IFischer et al.l |200 j) had values which were larger 
than those from most other sources by integer factors. This 
is due to the use in HITRAN of a different convention for 
the treatment of nuclear spin states. The convention in as- 
trophysics is not to count nuclear spin states as distinct 
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states in equation [S] As stated by llrwini l|l98lh . "the sta- 
tistical weights are divided by ... the product of the nuclear 
spin statistical weights". It is clear fro m the desc r iption s 
of the IflTRAN partiti on fu nctions iniFischer et al.1 (|2003h . 
ISimeckova et al.] l|2006l ) and lLaraia et ah 1 20 111 ) that a dif- 
ferent convention is adopted in HITRAN, and some other 
spectroscopy literature, where the corresponding statistical 
weights include nuclear spin states. 

The consequence is that on the HITRAN convention, 
partition functions Qhit are numerically larger than astro- 
physical partition functions Qast by a factor defined as fol- 
lows: 



Table 3. Continuum absorption processes in VSTAR 
Process Reference 



H~ free-free 
H~ bound-free 
H bound-free and free-free 
free-free 
Ha CIA 



-'2 

Ha 



Ha - He CIA 



Oa 
Ha 



Oa and Oa 
Na CIA 



Na CIA 



Bell fc Berrington (1987'! 

Wishart (1979) 

Gray (200^ ) ppl49-154 

Bell (1980) 

Borysow ( 2002) 

Borysow et ah (2001) 

Borysow et al. (1997, 1989) 

Boryso w &: Frommhold (_1989) 

Smi th fc Newnham (1999 . 2000) 

McKellar (1989) 



Qhit(T) = (3ast(r) ]^(2/j + 1) 



(9) 



where Ij are the nuclear spins of the n atoms contained 
in a molecule. With 7 = for carbon and oxygen, ^/a for 
hydrogen and 1 for nitrogen, these factors becone 4 for HaO, 
16 for CH4 and 24 for NH3. 

In many cases (e.g. equation [2ll partition functions are 
divided such that these integer factors will cancel out. How- 
ever, this is not the case when line intensities are calculated 
from Einsten A coefficients using equation [S] In this case 
it is important that the statistical weights (gs) and parti- 
tion function use the same convention. For example for HaO 
lines, if an astrophysical partition function is used the sta- 
tistical weights for the ortho and para states of water are 
/4 and ^/4, but if the HITRAN partition function is used 
they are 3 and 1. 



2.5.3 Atomic Partition Functions 

Partition functions for atoms and atomic ions in VSTAR are 
calculated using a modified version of the PFSAHA subrou- 
tine taken from the ATLAS9 software of Kurucfl 



2.6 Continuum Processses 

A number of continuum absorption processes are included 
in VSTAR as described in table [S] 



2.6.1 Bound-Free and Free-Free Continuum Absorption 

Bound-free and free-free absorptions of a number of species 
are important in stellar atmospheres, and some of these pro- 
cesses remain significant at the cooler temperatures of hot 
Jupiters and brown dwarfs. The processes included in VS- 
TAR are listed in table [S] These processes are sufficient for 
the range of temperatures and wavelengths currently stud- 
ied using VSTAR. A wider range of continumm processes 
would need to be incorporated to extend the applicability 
of VST AR to higher temperatures and into UV wa velengths 
(see e.g. ISharp fc Burrowsll2007i : ICustafsson et ahllioOSi ). 



2.6.2 Collision Induced Absorptions 

Collision induced absorption (CIA) due to Ha - Ha pairs 
is important in the solar system giant planets where it 
shows up a distinct strong spectral feature at 2.1 /im (see 
figure [S] for the spectrum of Jupiter in this region). This 
absorption is also important in brown dwarf, hot jupiters, 
and M dwarfs. We incorporate this absorption by interpo- 
lating in tables provide by Borysow[f|. Three different ta- 
bles are needed to cover t he fu ll temperature range using 
calculations bv iBorvsowl (l2002h (6 - 350 K and 400 - 
1000 K) and iBorvsow et al.1 (|200ll ) (1000 - 7000 K). Ha 
- He colli sion induced absorpt i on is also included using 
data f rom Bor vsow et al. I (ll997l ). lBorvsow et al.1 (|l989l ) and 
iBorvsow fc Frommholdl (|l989l ). 

Other collision induced absorptions included in VS- 
TAR are absorptio n in the Oa near infrared bands 

( Smi th fc NewnhamI 1 19991 . I2OO0I ) which is needed for the 
Earth atmosphere, and Ha - Na collision induced absorp- 
tion which occurs in the Titan atmos phere and for which 
data is available from iMcKeUail l| 19891 ). 



2.7 Scattering 

Scattering processes are often treated in stellar atmosphere 
codes as just an additional source of opacity which can be 
added to that due to gas absorption. However, VSTAR uses 
a more rigorous treatment that requires a full description of 
scattering. In addition to the optical depth contribution, we 
need to know the single scattering albedo, uj and the phase 
function, P{6), that describes the angular distribution of 
scattered light. 

Single scattering albedo is defined as: 



(10) 



where ae, a a and a a are the extinction, absorption and 
scattering cross sections per particle. 

2. 7. 1 Rayleigh Scattering from Molecules 

VSTAR includes models for the Rayleigh scattering due to 
air, H, Ha, He, Na and CO3 or any mixture of these gases. 
The wavelength dependence of the scattering cross section 



http: / /kurucz. harvard.edu/programs/atlas9 



http : / / www. astro. ku. dk / aborysow /programs 
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of H2 is from iDalgarno fc Williamsl J 19621) . Those of H and 
He are from Dalgarno as cited by Kurucz l|l970l ). 

For the other gases the scattering cross section (a) 
is derived from the refractive index of the gas using 
l|Hansen fc Travis! 1 1974 ): 



87r^(n(A)^ - 1) = 



Fk 



(11) 



where n(A) is the wavelength dependent refractive index, A 
is wavelength and N is the number density (molecules cm"'^) 
of the gas measured at the same temperature and pressure 
as the refractive index. Fk is known as the King factor and 
is given by: 



Fk = 



6-75 



(12) 



where 5 is the depolarization factor. For ai r we use the re- 
fractiv e index wavelength dependence f rom Peck fc Reedeij 
l| 19721 ) and the King fa ctor from Youn^ (|l98ll ). For CO^ the 
refracti ve index is from Old et al.l ( 197l|) and the King factor 
is from ISneep fc Ubachd (|2005D ." For N2 the refractive index 
is fromlCox (2000). 

Rayleigh scattering is pure scattering so it has a single 
scattering albedo of one, and the phase function is: 



P {9)^^(1 + cos' e) (13) 



2.7.2 Scattering from Particles 

Scattering from liquid or solid particles suspended in the at- 
mosphere (variously called clouds, aerosols or condensates) 
occurs in all planetary atmospheres in the solar system, and 
in brown dwarfs and M stars. 

Particle scattering can be modelled using Lorenz-Mie 
theory, which provides a solution to Maxwell's equations 
that describes scattering of light from a homogenous sphere 
with a complex refractive index n + ik, and a size param- 
eter X = 2'Kr/X where r is the radius of the sphere, and 
A is the w avelength. VSTAR uses a Lorenz-Mie scattering 
code from iMishchenko et all l|2002l ) that models scattering 
from a size distribution of spherical particles. The particle 
size distribution can be described by a number of functional 
forms including power law, log normal and gamma distrbu- 
tions. VSTAR can handle up to 20 particle modes, each with 
its own wavelength dependent refractive index, particle size 
distribution, and vertical distribution in the atmosphere. 

The scattering code generates the extinction and scat- 
tering cross sections per particle (cte and as) which give the 
single scattering albedo according to equation 1101 

It also provides the phase function P{6) both as an array 
of values for different scattering angles, and as an expansion 
in Legendre polynomials 



P{e) = q"P, (cos ^ 



(14) 



which is the form required for the discrete-ordinate radiative 
transfer code described in section \2l8\ 

If the particles are genuinally spherical (e.g. liquid 



droplets) then the phase function derived from the Lorenz- 
Mie code is appropriate. Sometimes however, spherical mod- 
els are used simply as a generic representation of solid parti- 
cles which probably have a range of irregular shapes In this 
case the full spherical particle phase function is probably 
not a good representation. It will include angular structure 
in the form of features such as rainbows and glories, that 
are specific to the particles' sphericity, and disappear for 
non-spherical particles (Bailey 2007) . In this case a better 
repre sentation may be the Heny ey-Greenstein phase func- 



repre sentation may be tne tleny i 
tion (|Henvev fc Greensteinlll94ll ') 



Phg{0) = 



1-9' 



47r(l + g^~ 2g cos 61)3/2 
where g is the asymmetry parameter defined by: 



g={cose)= cos 6'P(6') 27r sin ^ 
Jo 



(15) 



(16) 



Here g is between —1 and 1, and is positive for for- 
ward scattering and negative for backward scattering. The 
asymmetry parameter is also provided by the Lorenz-Mie 
scattering code. 

As well as using Lorenz-Mie theory, VSTAR can make 
use of precalculated scattering parameters derived using 
other methods. Scattering from non-sphe rical particles can 
be modelled usin g T-matrix methods (IWatermanl Il971 



|^^henk3|T993). Codes available from IMishchenko et aP 
l|2002t ) can be used to calculate the scattering properties of 
particles with spheroidal or cylindrical shapes or shapes de- 
scribed by Chebyshev polynomials in either random or spe- 
cific orientations. However, T-matrix techniques are gener- 
ally limited to small size parameters and can be very slow to 
compute. Some exampl e s of a pplic ations of thes e techn iques 
can be found in lBailevI l|2007l ) and lBailev et al.1 (|2008l ). 

The aerosols that make up Titan's stratospheric haze 
have optical properties that can be modelled by fractal ag- 
gregrates compos ed of man y spherical particles. Models de- 
scribed by iTomasko et al.l |2008i' ) give the scattering prop- 
erties of the Titan aerosols including optical depth, single 
scattering albedo and phase function in a form suitable for 
use with VSTAR. 



2.8 Radiative Transfer 

The data on absorption and scattering derived by the meth- 
ods described in sections 12.21 to 12.71 are finally combined to 
provide the inputs needed for solving the radiative transfer 
equation. VSTAR is structured in a modular way so that 
a number of different approaches to radiative transfer can 
potentially be used. However, for most of our work we have 
used the DISORT package ((stamnes et al.lll988l ). DISORT 
is a robust computer implementation of the discrete ordi- 
nate method for radi ative transfer originally developed by 
IChandrasekhEul l|l960l) . 

DISORT models the radiative transfer in a multiple 
layer medium, including the processes of absorption, multi- 
ple scattering and thermal emission, and allows a reflecting 
surface at the lower boundary and a direct beam source (e.g. 
the Sun) illuminating the top of the atmosphere. 
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The radiative transfer equation solved by DISORT has 
the form: 



' dr 



(17) 



where is the monochromatic radiance (sometimes referred 
to as intensity or specific intensity) at frequency v, and is a 
function of optical depth r, and direction fi, </>, where fi is 
the cosine of the zenith angle, and 4> is the azimuthal angle. 
The source function Sv is given by: 



471- JO J^l 

+ il-uj)B,{T) 

+ ^;^-P(m, 0; MO, (t>o) exp {-r/fio) 



(18) 



where the first term describes scattering of radiation into 
the beam from other directions according to single scattering 
albedo zu and phase function P(/i, </>; fi' , (j>'), the second term 
is thermal emission, with Bv(T) being the Plank function 
and the third term is direct illumination of the atmosphere 
by an external source with flux ^oFy and direction /io,<^o 
(e.g. the Sun). 

The inputs required by DISORT are the vertical optical 
depth (Ar), single scattering albedo (ro) and phase function 
moments (a", the coefficients of the Legendre polynomial 
expansion of the phase function as defined in equation I14p 
for each atmospheric layer. The temperature at each level is 
also needed for thermal emission calculations. 

The values of these are obtained by summing the con- 
tributions of all relevant processes for each layer of the at- 
mosphere at each wavelength. The combined absorption due 
to all molecular and atomic line processes as well as related 
continuum processes (sections 12.21 to 12. 6|) is called Arjas. 
Rayleigh scattering (section I2.7.1|) contributes an optical 
depth Array which is pure scattering optical depth. Each 
of the p modes (p = 1 to Pmax) of scattering particles (sec- 
tion [2I7]21) provides both a scattering optical depth At^^^^^ 
and an absorption optical depth At^^^. 

The total optical depth for a layer is: 



Ar = Ar,a, + Array + ^ (Arf,,,, + Arf.J (19) 

The combined single scattering albedo for a layer is the 
scattering optical depth divided by the total optical depth 
so is given by: 



(20) 



and the combined phase function moments (a^ for s = Q 
to Smax) are the phase function moments for the individ- 
ual scattering components, the Rayleigh scattering (aray) 
and the p particle modes {a"^), weighted according to their 
contribution to the scattering optical depth. 



'-'ray'-^ ' ray ^ Z^p—l " '-^ ' acatt 
Array + X]p=l'' ^"^scatt 



(21) 



DISORT approximates the angular distribution of the 
radiation field, by replacing the integral term in the source 
function with a sum over a number of discrete zenith angles 
(or streams) according to the Gaussian quadrature rule. The 
number of streams determines how accurately the angular 
distribution of radiance is calculated. For many purposes 
where the variation with angle is reasonably smooth, we have 
found 8 streams (four upward and four downward) to be ad- 
equate. However, the number of streams can be increased as 
necesssary to provide more accurate representation of angu- 
lar structure, at the cost of increased computation time. 

DISORT provides a number of outputs. In astronomical 
cases we are mostly interested in the upward radiance and 
its angular dependence at the top of the atmosphere. The 
radiance can be determined for any emission angle, and for 
any azimuth relative to the illuminating source. It is also 
possible to determine the radiation field at the surface of 
the planet or at intermediate levels in the atmosphere. 

The radiative transfer equation only holds for 
monochromatic radiance. Thus to obtain the spectrum of a 
planet the radiative transfer calculation must be repeated for 
each wavelength, and the wavelength steps made sufficiently 
small to resolve the spectral lines. We normally use a set of 
points that are equally spaced in wavenumber, and with a 
spacing chosen to be just sufficient to resolve the narrowest 
absorption lines encountered in any layer of the model. Typ- 
ically this requires several hundred thousand spectral points 
to be calculated. 



2.9 VSTAR Operation 

We can now outline the full procedure for using VSTAR to 
calculate the spectrum of a planet, brown dwarf or star. The 
first steps in the process are different depending on whether 
we are modelling a solar-system planet (including the Earth) 
or an exoplanet, brown dwarf or star. 



2.9.1 First Step - Solar System case 

In the solar-system planet case we start with a measured 
profile for the temperature and pressure as a function of 
height. This is available for all planets in the solar system, 
either from radio occultation observations, or from entry 
probes or other in-situ data. For Earth there are a number 
of "standard atmosphere" profiles available and several of 
these are built-in to VSTAR. Some other planets have stan- 
dard atmospheres, e.g. the Venus I nternational Reference 
Atmosphere (VIRA. [Seiff et al.lll985D 

The mixing ratios of the various gases and their dis- 
tribution with altitude are also generally reasonably well 
known from in-situ measurement or past spectroscopic anal- 
ysis. Thus for any such planet we can start with an atmo- 
spheric profile which consists of a number of layers (typi- 
cally 30 to 60) with the pressure, temperature, and mixing 
ratios of absorbing gases specified for each layer. If clouds or 
aerosols are included in the model their vertical distribution 
is specified in terms of the aerosol optical depth at a specific 
reference wavelength, for each particle mode at each layer. 
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2.9.2 First Step - Giant Exoplanet, Brown Dwarf or Star 

In these cases we don't have the direct measurements of 
profile and composition available for solar system plan- 
ets. Once again we start from a pressure-temperature pro- 
file for the atmosphere. The other input required for the 
model is the elemental composition of the atmosphere. 
For this we normally assume either a solar composition 



fe.g.lGrevesse. Asplund fc Sauvai2007l : lAsplund et al.ll2009l : 
iLodders. Palme fc Gailll2009l ) or a modified solar composi- 
tion for a different metallicity. 

Then for each layer of the atmosphere we use the ICE 
chemical equilibrium model (section 12. ip to determine the 
equilibrium composition of the layer in terms of molecules, 
atoms and ions. ICE also predicts what condensates are pro- 
duced at different levels in the atmosphere, and this can be 
used as a guide to the addition of clouds to the model which 
are specified in the same way as in the solar system case. 

Note that VSTAR is not an atmospheric structure 
model, and so cannot be used to determine structures that 
are self consistent with the radiative transfer (although this 
capability may be added in the future). Currently, therefore, 
we normally work from a pressure-temperature structure 
determined from another stellar or exoplanet atmosphere 
model and VSTAR is used as a spectral synthesis model. 



2.9.3 Second Step - Computing Layer Absorption and 
Scattering Properties 

In either case we now have a specification of the atmospheric 
layers with their pressures, temperatures, chemical compo- 
sition and aerosol optical depth. We now calculate for each 
layer at each modelled wavelength (typically there are sev- 
eral hundred thousand wavelength points) the absorption 
and scattering properties. For line absorption this involves 
calculating the line profiles of all the spectral lines of all rel- 
evant species that fall within the wavelength range of the 
model (including lines outside the range where the far wing 
contribution may be significant) and adding their contribu- 
tions into the gas optical depth of the layer. 

For particle scattering (clouds and aerosols - section 
I2.7.2P there are two options. Lorenz-Mie theory can be used 
directly to calculate the scattering properties of each parti- 
cle mode at each layer for a set of wavelengths. These values 
are then spline interpolated to provide data for each individ- 
ual wavelength point. Alternatively precalculated scattering 
properties can be used and interpolated in the same way. 




2,2 2.3 
Wavelength (jum) 



Figure 3. Comparison of Earth atmosphere zenith transmission 
calculated with VSTAR and RFM in the 2-2.5 /im window. Lower 
panel — Transmission calculated with VSTAR. Upper panel — 
Percentage difference (VSTAR - RFM). 



looking at reflected solar light. For a star we generally want 
the flux spectrum (the radiance integrated over all angles) 
which is equivalent to what would be seen from an unre- 
solved star at a large distance. However, the angular depen- 
dence of the radiance also allows the investigation of limb 
darkening of a star. VSTAR can also calculate the trans- 
mission spectrum from the top of the atmosphere to the 
surface. This is used in particular for the Earth atmosphere 
case where such data can be used to model tellu ric correc- 
tions to astronomical spectra l|Bailev et al.ll2007i l. 



TESTS OF VSTAR 
OTHER MODELS 



COMPARISON WITH 



We have tested VSTAR by comparing its predictions with 
a number of other models. 



2.9.4 Third Step - Radiative Transfer Solution 

The next step is to combine the absoprtion and scattering 
properties of each layer as described in section [2]8] to provide 
the inputs needed for DISORT, and perform the radiative 
transfer solution for each wavelength point. At this stage 
we specify the boundary conditions, which can include a 
reflecting surface at the base of the atmosphere, and an il- 
luminating source (the Sun or another star) at a specified 
zenith angle. 

VSTAR can calculate a number of different types of 
spectra. In the case of a solar system planet we are ususally 
interested in the spectrum of the radiance at a specific point 
on the planet's disk, or the radiance factor (I/F) if we are 



3.1 Comparison with RFM (Reference Forward 
Model) 

The Reference Forward Model (RFMi3 

is a line- by- line 

model for the Earth atmosphere developed at Oxford Uni- 
versity for the MIPAS instrument project on the Envisat 
satellit e. RFM is itself a development of the earlier GENLN2 
model l|Edwardslll993 ). RFM has been included in a num- 
ber of intercom parison studies of Earth atmosphere radiati ve 
transfer codes (|Tiemkes et al.ll20o3 : ISaunders et aLlbOOTl ). 



* http://www.atm.ox.ac.uk/RFM/ 
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We compared VSTAR with RFM (version 4.25) by us- 
ing both codes to calculate the Earth atmosphere transmis- 
sion spectrum at the zenith in the 2 to 2.5 wavelength re- 
gion. The atmospheric profile was the standard mid-latitude 
summer atmosphere taken from the RFM website. Absorp- 
tion due to H2O, CO2, O3, N2O, CO, CH4 and O2 were 
included in the models with line data taken from the HI- 
TRAN 2000 database (the version of HITRAN supported 
by RFM). 

The transmission spectrum from VSTAR and the per- 
centage difference between transmission calculated by RFM 
and VSTAR are shown in figure[3] The mean difference (VS- 
TAR - RFM) in transmission is -0.05%, and the RMS dif- 
ference is 0.09%. There is no significant trend in the dif- 
ferences with wavelength. These results are basically a test 
of the line-by-line absorption calculations and show the two 
models are in good agreement. Residual differences between 
the two codes appear to be due to differences in the line 
shape model used, and to differences in the way layer proper- 
ties are interpolated from the levels listed in the atmospheric 
profile. 



3.2 Radiative Transfer Benchmarks 

We have used VSTAR to reproduce a number of 
benchmark probl e ms in radiative transfer described by 
iGarcia fc Siewer3 (|l985l ). The first problem considered is 
a uniform Mie scattering atmosphere containing spherical 
particles with size parameter 2 and index of refraction 1.33. 
The atmosphere has a total optical depth r = 1, and the 
single scattering albedo zo — 0.95. The atmosphere is illu- 
minated by a source at a zenith angle of 60 degrees and 
has a non-reflecting surface at its base. The problem was 
implemented using VSTAR's own Mie scattering code to 
calculate the phase function, and the results for the top 
of atmosphere upward radiances and bottom of atmosphere 
downward radiances at an azimuth of zero are given in table 
131 and compared with th e results summed from t he fourier 
series coefficients given in iGarcia fc SiewertI l| 19851 ). The VS- 
TAR res ults, calculated wi t h 16 s treams in DISORT, agree 
with the IGarcia fc SiewertI 11985") results to five significant 
figures in all cases, and at over half the zenith angles agree 
to within ±1 in the sixth figure. Similar agreement is found 
for results at other azimuths and optical depths. 

Another problem we have calculated with VSTAR is 
one of the test problems posed by the Radiation Commis- 
sion of the Internatio nal Associatio n of Meteorology and 
Atmospheric Physics l|Lenoblel ll977l V This involves an at- 
mosphere with T = 1, zu = 0.9 with a "Haze L" size distri- 
bution of scattering particles as defined bv iLenoblel (|l97'if ). 
The atmosphere is illuminated by a source at a zenith an- 
gle of 60 degrees. This problem has a more complex phase 
function , and with 16 DISORT s treams the benchmark re- 
sults in IGarcia fc SiewertI (Il985l ) are only matched to an 
accuracy of ~1 per cent. However, increasing the number 
of DISORT streams to 64 p roduces results that agree with 
the IGarcia fc SiewertI (|l985t ) results to ±1 in the sixth sig- 
nificant figure for the upward radiances at the top of the 
atmosphere. 



Table 4. Radiances for Mie scat tering benchmark problem com- 
pared with values from iGarcia & : Sicwcrt (1985). Positive values 
of fj, are upward radiances at the top of the atmosphere and neg- 
ative values are downward radiances at the bottom of the atmo- 
sphere, both for an azimuth of zero relative to the illuminating 
source. 



H = coszd VSTAR Garcia &i Siewert Difference 
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3.3 Comparison with MARCS 

The MARCS stellar atmosphere code l|Gustafsson et ahl 
I2OO8I ) is a good choice for a comparison with VSTAR, since 
it uses modern abundances and opacities, similar to those 
used in VSTAR, and has been used to generate a large grid 
of model atmospheres with detailed information available 
on the model structure. The comparison shown here is for 
a plane parallel model with Teg = 2500 K, g = 1000 m s~^ 
(logg = 5 in c.g.s units), and solar metallicity. To make the 
comparison we take the temperature as a function of gas 
pressure from the MARCS model, and use this as the input 
for calculating an equivalent VSTAR model. 

Molecular species included in the VSTAR model were: 
H2O (BT2), CO, CaH, MgH, FeH, CrH, Ti O and VO. Lines 
of al k ali metals were ta ken from VALD l|Piskunov et al.l 
1 19951 : iKupkaet al.] 1 19991) . and other atomic species were 
from the Kurucz line lists. 

Figure m shows some comparisons of data from MARCS 
and VSTAR. These include the height above the lowest level 
modelled (determined from hydrostatic balance), and the 
mixing ratios of several important species. On these plots 
crosses are the values from MARCS and the lines are the 
values from VSTAR. In most cases the agreement is excel- 
lence. The behaviour for TiO is however quite different in the 
cooler layers. This is because the VSTAR chemical model in- 
cludes condensate formation, whereas MARCS only includes 
gas phase chemistry. Small differences in the CO and H2O 
mixing ratios arise from the same cause. The difference in 
TiO does not have much effect on the strength of the TiO 
bands in the spectrum as these are mostly formed at deeper 
levels. 

The flux spectrum from VSTAR is compared with that 
from MARCS in the bottom right panel of flgure |4l The 
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Figure 4. Comparison of MARCS and VSTAR for a stellar atmosphere model with T^g = 2500 K, logg = 5 and solar metallicity. The 
height, temperature and mixing ratios of several important species are plotted as a function of gas pressure. On these plots crosses are 
the values from MARCS, and the line shows the values from VSTAR (only the temperature profile is taken from the MARCS model. All 
other quantities are calculated independently by VSTAR and MARCS). The fiux spectra from MARCS and VSTAR are shown in the 
bottom right panel, with the MARCS data offset by 1.5 xlO®. 



fluxes provided with MARCS are described as "rough esti- 
mates of the surface fluxes" and "are NOT synthetic spec- 
tra" . Nevertheless it is clear that the values of the fluxes are 
very similar and the same absorption fetaures are present 
in both spectra at similar depths. A comparison of similar 
models with observations of an M dwarf spectrum are pre- 
sented later. 



4 COMPARISON WITH OBSERVATIONS 

In this section we present a brief comparison of VSTAR 
models with observed spectra for a range of objects. Full 
details on these and other results will be given elsewhere, 
but they are presented here to show the wide range of as- 
tronomical objects that can be successfully modelled us- 
ing VSTAR. The spectra of solar system objects are from 
our own observations with the Anglo- Australian 3.9m tele- 
scope at Siding Spring o bservatory, and its IRIS2 instru- 
ment (|Tinnev et aLll2004) . and from the NIFS instrument 
l|McGregor et al.ll2003l ) on the Gemini North 8m telescope 



at Mauna Kea Hawaii. The brown dwarf and stellar spec- 
tra are taken from the NASA Infrared Telescope Facility 
(IRT F) spectral hbrary (|Cushing et al.ll2005l : iRavner et al] 
I2OO9I ) and are taken with the SpeX instrument on IRTF at 
a resolving power of R'-^2000. 

The spectra considered here are all in the near-infrared 
spectral region (1 - 2.5 fim). For late type dwarfs this is 
where the flux peaks, and for spectroscopy of planets, this is 
a region containing many interesting rovibrational molecular 
bands. While this is the wavelength region VSTAR has nor- 
mally been applied to, there is nothing inherent in VSTAR 
that restricts it to this region. VSTAR should, in principle, 
be usable from UV to microwave wavelengths. 



4.1 Solar System Objects 

Figure [5] shows a number of comparisons of observed spec- 
tra with VSTAR models for solar system atmospheres. 
The Earth atmosphere transmissi on comparison (top left 
panel) is des c ribed in more detail in iBailev et all l|200'iD and 
iBailev et all (|2008l ) and uses the observed spectrum at a re- 
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Figure 5. Comparison of observed spectrum of solar system objects with VSTAR model spectra. At top left is shown the spectrum of a 
G2V star seen through the Earth's atmosphere in the region of the O2 a-X absorption band, compared with a VSTAR model for a solar 
spectrum with Earth atmospheric transmission applied. The other three panels show observed reflectance spectra of Titan, Jupiter and 
Uranus compared with VSTAR models. The observed spectrum in each case is offset upwards relative to the model spectrum. 



solving power of R A/AA) ~ 2400 of the G5V star BS 996 
with the 3.9m Anglo-Australian Telescope (AAT) and its 
IRIS2 instrument. This spectrum is compared with a model 
of a solar spectrum as seen through a modelled Earth at- 
mosphere transmission spectrum calculated using VSTAR. 
Of the spectral features only that near 1.28 /im is a stellar 
line. Other features are atmospheric absorptions of O2, CO2 
and H2O, with the strongest feature being the O2 a-X band 
at 1.27 lira and its associated collision induced absorption. 
The data and model agree to better than 1%. 

The bottom left panel shows the reflected light spec- 
trum of Jupiter in the near-infrared K band (2.04 - 2.38 
ixm) also observed with the AAT and IRIS2. The VS- 
TAR model used here is described in more detail in 



iKedziora-Chudczer fc BailevI |20l3). The absoprtion fea- 
tures present in this spectral region are due to methane, 
and to the collision induced absoprtion of H2-H2 (see sec- 
tion I2.6.2p . with the latter producing the broad smooth 
feature centred on 2.11 /xm. Stratospheric and tropo- 
spheric clouds are included in the model with the de- 
tails of the cloud properties and optical depths given in 
iKedziora-Chudczer fc BailevI (|201ll '). 

The right hand panels of figure [S] show spectra of Ti- 
tan and Uranus in the spectral region covering the 1.55 /xm 
methane "window" that lies between strong methane abso- 
prtion bands. Modelling of the spectra of these objects at 
this wavelength using line-by-line methods has only become 
possible very recently with the availability of the latest labo- 
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Figure 6. Comparison of observed spectrum of the T4.5 dwarf 
2MASS J055591915-1404489 witii a model spectrum calculated 
witli VSTAR for T^g = HOOK and log 3 = 5. The observed spec- 
trum is offset from the model upwards by 1 (x 10^^^) 



Figure 7. Expanded view of two sections of the spectrum shown 
in figure |6] The observed spectra are offset from the model up- 
wards by 0.5 (x 10"") 



ratory spectral line data for me thane as des c ribed in section 
in particular the data of I Wang etaH (120111 ). The Ti- 
tan spectrum was obtained with the Near Infrared Integral 
Field Spectrometer (NIFS) on the Gemini North 8m tele- 
scope, a,nd the data and VSTAR model ar e described more 
fully bv lBailev, Ahlsved fc Meadows! (120111 1. The absorption 
lines in this region are mostly weak lines of CII4, but there 
are also lines of CH3D and CO. A model for Titan's aerosols 
based on that of Tomasko et al. I dioos i is also used in this 
analysis. The models enabled Titan's D/H r atio and CO 
abund ance to be determined as described in iBailev et al.l 
(|201ll ). 

The Uranus comparison in the lower right panel uses 
data from IRIS2 on the AAT. The VSTAR model spectrum 
is calculated from the new low temperature methane line 
data as described in section [2.2.6l The model includes clouds 
with a mean particle size of 1 nm, and with the main cloud 
layers being at 2 - 3 bars and at 6 - 10 bars, which is similar 
to clo ud d istributions derived by ISromovskv. Irwin fc Frvl 
l|2006t ) and llrwin. Teanbv fc Da^ (120071 ). 

While not shown here VSTAR has also been used suc- 
cessfully to model the spectra of other solar system planets. 
VSTAR models fo r the V e nus n ight s ide are used in the an al- 
yses described bv lBailevI (|2009l ) and lCotton et all l|201ll ) 



4.2 Brown Dwarf 

Figure [6] shows an observed spectrum of the T4.5 
dwarf 2MASS J05591 915-1404489 taken from the 
IRTF spectral library (jCushing. Ravner fc Vaccal I2OO5I : 



I Ravner. Cushing fc Vaccal 20091 ) compared with a VSTAR 
model spectrum. The VSTAR spectra used for this compar- 



ison are based on brow n dwarf model P-T st ructures taken 
from Figures 8 and 9 of iBurrows et al.1 (|2006l ). The VSTAR 
models included molecular absorption due to H2O (BT2), 
CO, CH4, CaH, MgH, CrH, FeH, TiH. Lines of alkali 
metals are includ ed from VALD l|Piskunov et al.l Il995l : 
iKupka et al.l Il999f ) . The far wings of very strong sodium 
and potassium lines in the visible have a significant effect 
on the shape of brown dwarf spectra in the 1 //m region. 
We used a x factor model for the far wing shapes of these 
lines with an exponentially decrease between 500 and 7500 
cm~^ from the line centre, with the parameters adjusted to 
provide a good match to the data. More physically based 
models for the far wing shapes o f alka li me tal lines are 
descri bed bv lBurrows fc Volobuve^ (|2003h and lAllard et al.l 
(|2003l ). 

A number of models were tried with different tempera- 
tures and gravities, and best agreement with the data was 
found for a model with T^s = 1100 K, g = 1000 m s"^ (logp 
= 5 in c.g.s. units) and solar metallicity. This is in reasonable 
agreement with other results for this obiec tlStephens et al 
(^2009.) find Tcff = 1200 K, logg = 4.5 while Idel Burgo et al 
(I2OO9I ) find Toff = 1002 K, log 5 = 4.9 using high resolution 
spectra. 

The model spectrum agrees well with the observed spec- 
trum except in the wavelength region from around 1.6 - 
1.7 fim. As discused in section [2.2.61 the methane line lists 
currently available do not include hot bands in this region 
and therefore underestimate the total absorption, and in- 
correctly model the observed structure. This problem exists 
with all models for T-dwarf spectra. There is also a dis- 
crepancy between the model and observations in the 2.1 
region where the model is slightly too low. 

In other wavelength regions, where methane absorption 
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Figure 8. Comparison of observed spectrum of the M8 dwarf 
VBIO witli a model spectrum calculated with VSTAR using a 
MARCS model structure for T^h = 2700K and logg = 5. The ob- 
served spectrum is offset from the model upwards by 3 (x 10^^'') 



Figure 9. Expanded view of two sections of the spectrum shown 
in figure |8] The observed spectra are offset from the model up- 
wards by 0.4 (x 10~13). The upper panel shows primarily absorp- 
tion in the E-A band of FeH. The lower panel shows absorption 
due to CO and H2O. 



4.3 M Dwarf 



is not significant, the agreement between model and obser- 
vations is excellent. What might, at first sight, appear to 
be noise on the spectra, is in fact spectral structure due to 
many absorption lines. This is apparent from figure [7| where 
two wavelength regions (1.1 - 1.27 /im and 1.4 - 1.58 /^m) 
are shown on an expanded scale. 

The model used here did not include dust absorption 
and scattering. This is acceptable for objects as late as T4.5. 
For earlier type brown dwarfs, and in particular the L-dwarf 
class, the effects of dust are important. While VSTAR is 
quite capable of modelling the radiative effects of dust in 
such systems, there are currently considerable diferences 
in the assumptions about dust properties used in different 
models. Some models use large (~100 /im) dust particles 
llBurrows et a l. 2006) while others use small (sub-/im) par- 
ticles (|Allard et al. 2001 ) . For the purposes of this paper 
we have avoided the extra complexity of dusty brown dwarf 
models. 

Our models also assume equilibrium chemistry. Some 
brown dwarf models include a non-equlibrium treatment of 
the chemistry of some species such as CH 4, CO and H2O 
ISaumon et al.ll2003l : [Stephens et~aLll2009l ') which allows for 
the effects of vertical mixing in the atmosphere. While this 
could be implemented as a modification to our chemical 
model we have not done so for the results presented here. 
The effects on the near-I R spectrum appear to be small 
(|Hubenv fc Burrows! [2OO7I ') with more significant effects at 
4-14 /im wavelengths. 



Figure [8] shows an observed spectrum of the MS dwarf 
VB 1 (Gl 752B) taken from the IRTF spectral library 
(1 Gushing. Ravner fc Vaccall2005l : [r avner. Gushing fc Vaccal 
l2009h compared with a model spectrum calculated with VS- 
TAR. The model structur e was taken from the MARGS 
grid of model atmospheres (iGus tafsson et aDl2008 h and is a 
model for T^s = 2700 K, g = 1000 m s"^ (log 5 = 5 in c.g.s. 
units) and solar metallicity. In other respects the model is 
similar to that described in section 13.31 A number of dif- 
ferent effective temperatures were tested with the 2700 K 
model providing the best match to the observed spectrum. 

Figure [9] shows expanded views of two regions of the 
spectrum. The upper panel shows the 1.58 - 1.78 /xm re- 
gion. This is dominated by absorption lines in the E-A band 
of FeH, which are included in our model using the line list 
of iHargreaves et al.l (|2010l ). This band is essential to get a 
good model of the spectrum in this region. However, in- 
cluding this line list results in a small step in the modelled 
spectrum at 1.582 ^m (visible in figure|8} which is the short- 
est wavelength included in this list. This step is not seen in 
the observed spectrum, and suggests that there is more FeH 
absorption at shorter wavelengths that should be included. 
The lower panel of figure [9] is 2.25 - 2.5 pm region showing 
that the model does a good job of representing the detailed 
structure of the GO and H2O bands in this region. As for 
the brown dwarf case we have used a dust free model. It 
is quite possible that dust is present in such a late type M 
dwarf, and its inclusion might change our conclusions about 
the best fitting effective temperature. 
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The upper limit on the layer temperature for modelling 
stellar atmospheres using VSTAR is currently 6000 K set 
by the use of data from the JANAF tables in our chemi- 
cal model. Since a model usually requires a range of layer 
temperatures extending to about a factor of two above and 
below the effective temperature of the star being modelled, 
the maximum effective temperature is around 3000 K. 



5 CONCLUSIONS 

We have desribed the techniques used in the VSTAR code 
to calculate model spectra for objects including solar-system 
planets, brown dwarfs and M dwarfs. While there are other 
codes used for modelling the spectra of these objects, VS- 
TAR has a number of unique features. 

• It is capable of being used for a very wide range of ob- 
jects, ranging from the coolest solar system planets (Uranus 
with layer temepratures down to 60 K) to stars with layer 
temperatures up to 6000 K. 

• It has been tested by comparison with both the Earth 
atmosphere radiative transfer code (RFM) and a stellar at- 
mosphere code (MARCS). 

• Not only is a rigorous approach to radiative transfer 
used, but this has been tested against radiative transfer 
benchmarks and shown to be accurate to levels of 10~^ to 
10'^ 

An obvious application of VSTAR is to the modelling of 
exoplanet spectra. VSTAR has the proven ability to model 
both cool planets in our own solar system, and hotter ob- 
jects such as brown dwarfs and M dwarfs, and thus covers 
the full temperature range expected to be encountered in 
exoplanets. The full treatment of the angular dependence of 
scattering incorporated in VSTAR means that it can be used 
to predict the phase dependence of an exoplanet spectrum 
around its orbital cycle. 

Work is currently in progress on extending VSTAR to 
include polarized radiative transfer. The phase variation of 
polarization can potentially provide importan t information 
about an exoplanet 's atmosph ere or surface (|Bailevl [200?! : 
IStamll2008l : IZugger et al1l2010l l 
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